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I. INTRODUCTION 



Analysis of a three flavor neutrino oscillation fit to recent Super-Kamiokande data 

Christoph Meier* and Tommy Ohlsson^ 
Division of Mathematical Physics, Theoretical Physics, Department of Physics, Royal Institute of Technology, 

SE-100 44 Stockholm, Sweden 
(February 1, 2008) 

We have analyzed the most recent available Super-Kamiokande data in a three flavor neutrino 
oscillation model. We have here neglected possible matter effects and performed a fit to atmospheric 
and solar Super-Kamiokande data. We have investigated a large parameter range, where the mixing 
angles were restricted to < 9i < n/2, i = 1,2, 3, and the mass squared differences were taken to be 
in the intervals 10"" eV^ < Am^ < 10"^ eV^ and 10"* eV^ < AM^ < 10 eV^. This yielded a best 
solution characterized by the parameter values 6i ~ 45°, 62 — 10°, S3 ~ 45°, Ani^ ~ 4.4 x 10"^" eV^, 
and AM^ ~ 1.01 x 10"^ eV^, which shows that the analyzed experimental data speak in favor of 
0^ ' a bimaximal mixing scenario with one of the mass squared differences in the "just-so" domain and 

On , the other one in the range capable of providing a solution to the atmospheric neutrino problem. 

PACS number(s): 14.60.Pq, 14.60.Lm, 12.15.Hh, 96.40.Tv 
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-vj . Most of the analyses concerning neutrino oscillations have so far been done within the framework of the theory of 

^ ■ two flavor neutrino oscillation. In these scenarios, the observed deficit of solar electron neutrinos is usually explained 
f^ ] by means of i/e '^ i^fj. oscillations, whereas the lack of atmospheric muon neutrinos is interpreted as a consequence 
t^^ ■ of i*^ <-> Vr oscillations. The main advantage of the two flavor oscillation scenario is that the oscillation probability 
Cn ^ depends on only two parameters, the mass squared difference Am^ and the mixing angle 9. 

Recently, it was pointed out by several authors, that in some cases the interpretation of the experimental data in 
terms of two flavor oscillations might yield misleading and sometimes even wrong results (see e.g. Ref. M) and that 
a three flavor oscillation description is to be favored. Unfortunately, the oscillation probabilities depend in this case 
on five parameters, two mass squared differences and three mixing angles, which makes it quite hard to deal with the 
dependence of the probability functions on these parameters. 

The most popular three flavor neutrino oscillation scenario is the so called bimaximal mixing scenario H , where two 

O I of the mixing angles are maximal, i.e., they have values of about 45°, whereas the remaining one is restricted by the 

^ CHOOZ experiment to be rather small at least for certain parameter ranges H. The three flavor oscillation scenario 

can then be shown to decouple into two distinct scenarios involving two flavors, in which case it is of course much 

easier to perform a fit to the experimental data. Many authors B use the CHOOZ result to simplify the oscillation 

probability formulas for three flavors and then perform a fit of the two decoupled two flavor neutrino oscillation 

^ . , scenarios to the experimental data. 

j^ ■ In this paper, we will present a numerical fit within a three flavor neutrino mixing scenario to the zenith angle 

distribution of 850 day atmospheric neutrino Super-Kamiokande data and 708 day solar neutrino Super-Kamiokande 

data. However, we will not make any assumptions about the parameters on which the oscillation probability formulas 

depend, i.e., we fit the scenario in the most general case and for a large parameter range to the experimental data. 

The paper is or ganized as follows. In Sec. O, we summarize the basic features of three flavor neutrino oscillation 



theory. In Sec. Ill, we present the choice of experimental data and in Sec. IV, we discuss the minimization procedure. 
In Sec. M, we give the obtained solutions of the minimization p roble m, which are interpreted in Sec. Wl. The fits 



corresp ondin g to the solutions are discussed and compared in Sec. VH and then these solutions are tested for stability 
in Sec. |VIII|. Finally, Sec. IX contains a summary and our conclusions. 
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II. THEORY OF THREE FLAVOR NEUTRINO OSCILLATIONS 

In three flavor neutrino oscillation theory one assumes the neutrino states with definite flavor \i'a), en = e,^,T, to 
be linear superpositions of states with definite mass \i'i), i = 1, 2, 3, 

3 

|i/q) =^ [/;;>,), a = e,M,r. (1) 

i=l 

The unitary mixing matrix U is called Cabibbo-Kobayashi-Maskawa (CKM) matrix and can generally be parameter- 
ized by three mixing angles 9i, where i = 1, 2, 3, and three CP-violating phases. One of the latter were recently claimed 
to be observable by next generation neutrino detectors H and could lead to physically very interesting consequences. 
The other two can be shown to cause no physical effects. However, for this analysis these phases are too small to be 
important and we will neglect them in what follows. This implies Uai = U^^ and therefore we can write the CKM 
matrix in its (real) standard parameterization Iq] as 

C2C3 S3C2 S2 

U == I -.S3C1 - S1S2C3 C1C3 - S1S2S3 S1C2 I , (2) 

S1S3 - S2C1C3 -S1C3 - S2S3C1 CiC2^ 

where Si = sin 6i and Ci = cos 9i . The probability for detection of a neutrino of flavor /3 in a beam of neutrinos 
consisting exclusively of flavor a at the source of the beam is for three flavors given by 

PaP = Sap- 4: 2_^2_^UaiUi3^UajUi3j sin I , tt, /3 = 6, /i, T, (3) 

i=l j=l \ J 

where Sap is Kronecker's delta, L is the distance from the source to the detector, E is the neutrino energy, and Am| 
is the difference between the squares of the masses corresponding to the mass eigenstates |fi) and {vj). Since 

Am2i + Am32 + Am^g = (4) 

the mass squared differences are not linearly independent, and it is convenient to choose 

Am^ = Am2i, (5) 

^32) 

from which follows that 



AM^ = Am^2, (6) 



Amgi = -Ami3 = Am^ + AAf^ (7) 

So far we have considered the neutrino states to be plane waves, i.e., to have a definite momentum. Since the 
neutrino state is neither produced nor detected with a definite momentum or propagation length, one has to average 
over L/E as well as other uncertainties in production and detection. We will here follow closely Ref. H] and assume 
that these uncertainties to be described by a Gaussian average 



(Pep) = / Pap{x)fix)dx, (8) 

J — oo 

where 

and X = L/E. Inserting formula (3) for Pap yields 

3 3 

(Pap) ^Sap~2Y,Y. Ua^Up^UaJUpJ [l - cos(2/Am2,-) exp{-272( Am^^- )2}] , (10) 

i<3 



where I is given by 



and therefore related to the sensitivity of the experiment. (L) should here be inserted in meters, whereas (E) should 
be inserted in MeV. The parameter 7 is the so called damping factor and in accordance with Ref. [Q we will choose 

where AE and AL are the uncertainties in neutrino energy and propagation length, respectively. 
For large values of 7 

exp{-2-f'^{Amlf}^0 

and the oscillation term vanishes. The transition probability then becomes a constant dependent on the values of the 
mixing angles; the oscillation is said to become "washed out". To proceed, one now has to determine the values of 
the parameters L, E, AL, and AE for the experiments under consideration. 

III. CHOICE OF EXPERIMENTAL DATA 

We are going to consider three types of Super-Kamiokande data, multi-GeV e-like atmospheric neutrino data, 
multi-GeV /z-like atmospheric neutrino data, and solar neutrino data. 

A. Solar neutrino data 

The probabilities were taken from Ref. [Q for 14 data points in an energy range from 6 MeV to 12 MeV. The 
probability is almost constant for all 14 data points. The path length for solar neutrinos is given by the distance 
from the Sun to Earth {L ~ 1.44 x 10^^ m) and the uncertainty of the path length is here assumed to be negligible 
compared to L, i.e., 

- - 0. (13) 



The damping factor then becomes 



L AE 

7-^ — , (14) 



where AE = 0.5 MeV is the energy resolution of the experiment. 

B. Atmospheric neutrino data 

The atmospheric neutrino data are divided into atmospheric multi-GeV e-like events and atmospheric multi-GeV 
/i-like events. In both cases, the flux measurements are performed for five bins, i.e., for five different values 

cos 61 = -0.8,-0.4,0,0.4,0.8 

of the zenith angle 9, which is defined as the angle between the direction of the measurement and the axis through 
experiment and the middle of the Earth. Therefore, cos 9 = — 1 corresponds to straight upward going neutrinos, 
whereas for cos 9 = 1 the neutrinos are produced in the atmosphere right above the detector. 

The experimental data used are 850 day Super-Kamiokande data and Monte Carlo no-oscillation predictions for 40 
years generated lifetime for the five bins. The corresponding errors were taken from Ref. [g|. 

The results of the Super-Kamiokande experiment are usually presented by means of the "ratio of the ratios" 

7^ = i^(data)/^(AfC), (15) 



where $„ is the flux of neutrino flavor a. 

In Summer 1998, the Super-Kamiokande CoUaboration pubhshed experimental values for this "ratio of the ratios" 
of 7^ = 0.63 ± OM{stat) ± 0.05{sys) for sub-GeV events and 7^ = 0.65 ± 0.03(stat) ± 0.08(sys) for multi-GeV events 
y|, which was the first direct observation of neutrino oscillations. 

However, to be able to fit both, e-like and /i-like data to the measurements, we will use here the directly measured 
disappearance rates 



for e-like events and 



^- - ^:jmc) ^'^^ 

for /i-like events. The main disadvantage of this is the uncertainty in the expected neutrino flux, which was claimed 
to be up to 20% Q, due to the unknown flux of cosmic particles hitting the atmosphere. To handle this one can 
introduce an overall fiux uncertainty factor r] with values between 0.8 and 1.2. The actual neutrino flux for the neutrino 
flavor a is then given by $„ = rj^a- Thus, in absence of neutrino oscillations, the ratio between the theoretical and 
the measured fluxes would just be Ra = t], whereas otherwise one obtains 

Ra = , a = e,fi, T, (18) 

where Pafi is the transition probability from a to (3, see Sec. O. In what follows, we are going to neglect the r 
contribution due to too small production cross sections. The uncertainty in the fiux of cosmic rays is then an overall 
factor common to the ratios Ra, where a = e, fi, of the both remaining neutrino flavors. Denoting R = $^/$e, one 
obtains 



Re = ViPee + RP^.^.) (19) 



for the electron neutrino ratio and 



1 

'r' 



R^=^[ Pee + IjPmM ) (20) 



for the muon neutrino ratio. The value of R is theoretically well determined to be about i? ~ 3 at -E ~ 10 GeV |ll[| . 
The energy of the measured multi-GeV neutrinos is assumed to be E' = 10 GeV with negligible uncertainty compared 
to E and we set /S.E/E ~ 0. The path length for a neutrino that was detected in a bin corresponding to the zenith 
angle 9 is given by 



L{e) = \/r2 cos2 e + 2rd + d?- - r cos 6, (21) 

as is easily obtained from geometrical considerations. Here r is the radius of the Earth and d is the typical altitude 
of the neutrino production point in the atmosphere, which we assumed to he d = 10 km. The uncertainty in path 
length is mainly determined by Acosfl = 0.2 and therefore 



Ai = 
The damping factor 7 is thus a function of 9 



dL{9) 



9 cos 6* 



rL 

A cos 61 = A cos 9 (22) 

L + r cos 6* ^ ' 



f^ 



IV. THE MINIMIZATION PROCEDURE 

According to the last section, we have 24 data points, which we can use for the fit to the parameters of the theory: 
the five bins for e-like events, the five bins for /i-like events, and the 14 data points for the solar neutrinos. 
The function used to obtain the parameters is given by 



5 -, 5 



1 f 1 1 



1 






w 



[Ps««,«-Ps««,*(ea;p)]'L (24) 



where the first two sums run over the five bins and the third one over the 14 data points for the solar neutrinos. The 
weights have been taken as 



W 



Re,i{exp)/ARe,i{exp), (25) 

M:» = Rf_,^i{exp)/AR^^^{exp), (26) 

^s«n,i(ea;p)/APsM«.i(exp), (27) 

5 5 14 



Wtot = > W 



^^fi.i+^Vfsun.i- (28) 



i—1 i—1 i—1 

The mixing angles were constrained to the interval [0, 7r/2], whereas for the mass squared differences we assumed one 
of them to be in the range 10"" eV^ < Am^ < 10"^ eV^ and the other one in lO"'^ eV^ < AM < 10 eV^ Since 
these large parameter ranges yield no numerically stable solution, we minimized function (GJ) for 45 parameter ranges 
corresponding to all combinations of values for the mass squared differences between 

10"" eV^ < Ato^ < 10"^" eV^ 
10"^° eV^ < Ato^ < 10-9 eV^ 

10"^ eV^ < Ato^ < 10^2 eV^ 

for the small mass squared difference and 

10"'* eV^ < AM^ < 10"^ eV^ 
10"^ eV^ < AM^ < 10"^ eV^ 

1 eV^ < AM^ < 10 eV^ 

for the larger one. 

Following Ref. ||l|, we chose the following strategy to deal with the flux uncertainty factor ry: During the search for 



minima of function (24), we set it equal to unity, i.e., we neglected the flux uncert ainty . The solutions obtained in 



this way were then tested for their stability with respect to variations of 7], see Sec. VII] . 

To minimize the function x in Eq. (|2J) we generated A'^ = 10^ random values for each of the flve parameters in the 
corresponding range such that we obtained N values of the minimization function. After that we picked the parameter 
set that yielded the smallest value for the function and repeated this procedure n = 20 times. Each of the obtained n 
parameter sets then served as starting values for a deterministic minimization procedure, using a sequential quadratic 
programming method such that one obtains again n parameter sets as an end result. 

V. SOLUTIONS OF THE MINIMIZATION PROBLEM 

The solutions to the minimization problem can best be obtained by considering tables like Table E, where the 
best point values for the two mass squared differences for each of the investigated mass regions are shown. One can 
see, that for most of the mass regions at least one of the two mass squared differences was obtained on one of the 
region boundaries. In this case, we assumed that the value of the corresponding mass squared difference converges 
towards a better minimum in one of the neighbor regions. Thus, a given parameter region was only considered to 
contain a solution to the minimization problem, if both values of the mass squared differences were obtained within 
the boundaries of the region. Note, that the values obtained for AM^ are almost randomly distributed in the last 
column of Tabic |, i.e., in the range 1 eV^ < AAf ^ < 10 eV^. This is due to a too large value of the damping factor 7 
in Eq. (nw such that the minimization function becomes independent of this parameter. Only Attj^ and the mixing 
angles are in this case fitted to the experimental data. 



To keep overview over the various considered parameter regions one can consider Table ||. Here a region, where a 
value for AM^ was obtained within the boundaries, is marked with an "M" , whereas a region with a corresponding 
value for Am^ is marked with an "m" . One obtains in this way eight regions, marked "M m" in Table ||, where the 
mass squared differences are within the region boundaries. Seven of these regions are contained in the first two rows of 
Table |l| whereas the eighth one was obtained within the range 10"'' eV^ < Am^ < 10"^ eV^, 1 eV^ < AM^ < 10 eV^. 
The parameter values obtained from the minimization procedure with the mass squared differences constrained to 
these regions correspond to the possible solutions of the minimization problem. One notices that the regions marked 
with "Af ?7i" split Table || in three distinct areas, corresponding to three types of solutions, as we will see. 

The first of these interesting areas is the one containing the parameter regions in the upper left corner of Table 0, 
i.e., the nine regions in the range 10~^^ eV^ < Ato^ < 10~* eV^, 10"'' eV^ < AM^ < 10"^ eV^. In this range, there 
are all together five possible solutions, i.e., mass regions, where the mass squared differences were obtained within the 



region boundaries. To pick the best of the possible solutions for this type, consider Table [I], where the best point 
values for the minimization function (p3) arc shown. The smallest value of the minimization function was obtained 
in the region corresponding to lO^^*' cV^ < Ato^ < 10^^ cV^, 10^^^ cV^ < AM^ < 10~^ eV^. In all surrounding 
regions, we obtained values for the minimization function, which are larger, such that there seems indeed to be a 
minimum in the region specified above. Thus, the corresponding parameter set is a solution of the minimization 



problem, and we will denote it by Solution 1 in what follows. Considering the rest of Table |I|, one can see that this 



solution corresponds to the smallest value of the minimization function in the whole inve stiga ted parameter range, 



which means that it provides the best fit to the experimental data, as will be seen in Sec. VII. Note that the value 
of the minimization function in the parameter range containing Solution 1 is very close to the one obtained in the 
region 10"^" eV^ < Am^ < 10~^ eV^, lO"'' eV^ < AM^ < 10"^ eV^. Furthermore, the values for the mass squared 
differences obtained within these two regions are rather close to each other, and one can not really distinguish, if these 
two possible solutions correspond to the same minimum or to two distinct minima very close to each other. 

The second interesting area of Table contains six parameter regions in the range 10~^ eV^ < Am? < 10"'^ eV^, 
10^^ eV^ < AM^ < 10 eV^. Here a second type of solution is obtained in the region 10"'' eV^ < Am^ < lO^'^ eV^, 
1 eV^ < AAf^ < 10 eV^, i.e., in the last column, eighth row of Table |l|. We will denote this solution by Solution 2 in 



what follows. Considering now Table III, one can see that in one of the neighbor regions we obtained a smaller value for 
the minimization function. One has therefore only a "local" minimum inside the region containing Solution 2, and the 
values of the minimization function in the surrounding regions do not converge that obviously towards this region, as 
they did in the case of Solution 1. However, in all the surrounding regions at least one of the mass squared differences 
was obtained on one of the region boundaries, which is why we will consider the region 10~* eV^ < /S.m? < 10"'^ eV^, 
1 eV^ < AM^ < 10 eV^ to contain the second solution of the minimization problem. Comparing the corresponding 
value of the minimization function to the one of Solution 1, one sees that this second solution provides a less exact fit 
to the experimental data. 

Finally, the third interesting area is situated in the upper right corner of Table O, in the range 10^^^ eV^ < Arn^ < 
10~ eV , 10~ eV < AAf < 10 eV . Here two regions with both mass squared differences within the region 



boundaries are obtained. Consideration of Table HI tells us, that the smallest value for the minimization function 
was obtained in the parameter region 10"^" eV^ < Airi^ < 10~^ eV^, 1 eV^ < AAf^ < 10 eV^. Thus, we obtain a 
third type of solution of the minimization problem, denoted by Solution 3 in what follows. As in the case of Solution 
2, there are neighbor regions with smaller values for the minimization function, but they both have at least one of the 
mass squared differences on one of the region boundaries. Note that the value of the minimization function is about 
five times larger for Solution 3 than the one corresponding to Solution 1 , and roughly three times larger than the one 
corresponding to Solution 2. 

VI. INTERPRETATION OF THE OBTAINED SOLUTIONS 

The parameter values for the three solutions obtained from the analysis in the last section are shown in Table 



IV. Here the average value obtained from the n = 20 minimizations are depicted, together with the corresponding 
standard deviations as well as the best point values. Solution 1 actually corresponds to the most common three flavor 
neutrino oscillation scenario, with one of the mass squared differences in the "just-so" domain (lO^-'^eV^ — lO^^eV^) 
and the other one in the range of possible solutions of the atmospheric neutrino puzzle (lO^^eV^ — lO^^cV^). As 
for our solution, in this common scenario the values of the mixing angles are normally a set with a small value for 
02, and values around 45° for the other two mixing angles, 9i and ^3, i.e., one has bimaximal mixing. The three 
flavor oscillation scenario can in this case be shown to decouple into two independent oscillation scenarios involving 
two flavors Q. The two flavor oscillation scenario with the small mass squared difference is then commonly used 
to describe the solar electron neutrino deficit in terms of v^-v^^ oscillations, whereas the scenario with the larger 



mass squared difference is assumed to describe the atmospheric muon neutrino deficit in terms of v^-v^ oscillations. 
However, this approximation is only exact in the case O2 — > 0, and the nonzero value of this angle causes a mixing of 
these two decoupled two flavor scenarios. This influences in turn the values obtained for the mass squared differences, 
which might explain the deviation of our result for the large mass parameter AM^ = 1.01 x 10"'^ eV^ from the Super- 
Kamiokande result AAf ^ ~ 3.5 x 10^'^ eV^ H, which was obtained in the framework of a two flavor oscillation scenario. 
Note that we obtained this bimaximal solution without putting any restrictions on the parameters. Furthermore, the 
value we obtained for the second mixing angle 62 — 10° implies that sin^ 6*2 — 0.03, which is in accordance with the 
upper bound obtained from the CHOOZ experiment sin^ ^2 < 0.05 P|. A disadvantage of Solution 1 is the rather 



large standard deviation of the value obtained for Am^, see Table IV. 

Solution 2, on the other hand, is characterized by a value for the minimization function, which is twice as large as the 
one corresponding to Solution 1, but all parameters, except for the large mass squared difference, were obtained with 
very small standard deviations. As discussed earlier, the minimization function is in this mass range independent of 
AM^. This means, that one obtains random numbers for this parameter, which explains the large standard deviation 
of AM^. Only the mixing angles and Ato^ are in this case fitted to the experimental data. Here the small mass 
squared difference is lying in the range of possible solutions to the atmospheric neutrino puzzle, whereas the large 
one is compatible with the results of the Liquid Scintillator Neutrino Detector (LSND) experiment II4I, which yielded 
a value for the mass squared difference in the eV-range within the framework of a two flavor oscillation scenario. 
However, one can see from Table pA^ , that the value obtained for the second mixing angle is maximal, i.e., 82 — 45° 
and, unlike for the first solution, it is therefore in this case not possible to reduce the three flavor oscillation scenario 
to two oscillation scenarios involving two flavors. 

Finally, Solution 3 has the major disadvantage of a rather large value of the minimization function compared to 
the other two solutions, as discussed in the last section. The remarkably large standard deviation of the large mass 
squared difference can be explained as in the case of the Solution 2, but Solution 3 also shows a rather large standard 
deviation of the third mixing angle 9^. This solution is actually an "intersection" of the first two solutions, with one 
of the mass squared differences in the "just-so" range and the other one in the range obtained by LSND. All three 
mixing angles have values around 9i ~ 30°, where i = 1, 2, 3, which means that it is, like in the case of Solution 2, not 
possible to decouple the three flavor oscillation scenario into a pair of two flavor oscillation scenarios. 

VII. DISCUSSION OF THE OBTAINED FITS 

Figures |l| - 1| show the results of the fit. The best point probabilities obtained from the three solutions are depicted 
as well as compared to the corresponding experimental data including error bars. 

Let us first consider the results of the fit to the atmospheric electron neutrino data, shown in Fig. |l[ Here Solution 
1 and Solution 3 provide the best fits to the experimental data. Solution 2 shows larger deviations, especially for the 
two bins with the smallest zenith angles, i.e., the largest values of cosO. The neutrinos measured in these two bins 
come from right above the detector, which means that the background of cosmic particles is not shielded by the Earth. 
The experimental data corresponding to these zenith angles have accordingly the largest experimental deviations, and 
from Eqs. ( p5| ) - ( p7| ) one can see, that these bins therefore correspond to the smallest weights. This explains why the 
fit is less exact for these zenith angles. The fact that one obtains a zenith angle independent fit for Solution 3 can be 
understood from consideration of the oscillation lengths 



AttE 
|A?n 



L'^i = TTZIT-v «,i = 1,2,3, i^j, (29) 



y I 



where E is again the neutrino energy and the ^mfj are the mass squared differences. This solution has one of the mass 
squared differences in the "just-so" range and the other one in the eV-range. The oscillation length corresponding to 
the small mass squared difference is then far too long to make it possible for the experiment to see any variations in 
the oscillation probability, whereas in the case of the large mass squared difference the oscillations wash out before 



the neutrinos reach the detector, due to a too large value of the damping factor 7 in Eq. dlOf ). 

Considering next the fit to the atmospheric muon neutrino data, shown in Fig. H, one clearly sees that here Solution 
1 provides the best fit. Solution 2 has again larger deviations for the two bins with the smallest zenith angle, which as 
before can be explained by the fact that these bins have the largest experimental errors. Solution 3 provides again a 
constant fit to the experimental data, for the same reasons as the ones discussed above. It appears that this solution 
is definitely not capable of explaining the zenith angle dependence of the atmospheric muon neutrino flux measured 
by the Super-Kamiokande detector, and therefore this solution can be ruled out. 

Finally, all three solutions provide rather good fits to the solar neutrino data, see Fig. ||. Solution 1 and Solution 3 
are seen to yield a constant oscillation probability up to an energy of about 9 MeV. Above that energy the oscillation 



length becomes of the same order of magnitude as the Sun-Earth distance such that the oscillation probability 
obtained from the two solutions starts to vary with energy. Solution 2 yields a constant fit to the solar neutrino data, 
which can be understood, having in mind that this solution is characterized by one mass parameter in the range 
10""' eV^ < Ato^ < 10~^ eV^ and the other one in the eV-range. In both cases, the corresponding oscillation lengths 
are much shorter than the Sun-Earth distance, and all oscillatory effects wash out before the neutrinos reach the 
Earth. 

To summarize, only the Solution 1 provides a good fit to all three types of experimental data, whereas Solution 
3 can be definitely ruled out by the zenith angle dependence of the atmospheric muon neutrino data measured by 
Super-Kamiokande. It remains now to test the stability of the obtained solutions with respect to variation of the flux 



uncertainty factor introduced in Sec. [V 



VIII. STABILITY OF THE SOLUTIONS WITH RESPECT TO THE FLUX UNCERTAINTY 

As mentioned before, we performed the minimization procedure setting the flux uncertainty factor rj introduced in 
Eq. ( p^ ) and Eq. (gO|) equal to one. But this factor can vary between 0.8 < 77 < 1.2 and therefore one has to test the 
solutions obtained in the last section for stability with respect to variations of this flux uncertainty factor. 



We saw in the last section, that the solution we denoted by Solution 3 in Table IV can be ruled out by the zenith 
angle dependence of the atmospheric muon neutrino flux measured by Super-Kamiokande. We will therefore only test 
the stability of Solution 1 and Solution 2. To do this, we applied basically the same minimization procedure as the 



one described in Sec. IV, but this time for nine equidistant values of the flux uncertainty factor between 77 — 0.8 and 
rj = 1.2. Again, the angles were restricted to the intervals < 9i < 7r/2, where i — 1,2, 3. 

In the case of Solution 1, the small mass squared difference was constrained to 10"^" eV^ < Am^ < 10^^ eV^. To 
choose an interval for the large mass squared difference, one has to recall that in one of the neighboring regions of 
Solution 1 , we obtained parameter values, which are very close to those corresponding to this solution, as we mentioned 
in Sec. [V|. This could mean both that the two parameter sets correspond to two distinct minima or to the same 
one. We will here follow the latter assumption, i.e., we will assume that there is one minimum of the function (p3) 
somewhere between these two parameter regions. To avoid obtaining the mass squared differences on the boundary 
between the two regions, we widened the interval for AM^ to 10~^ eV^ < AM^ < 10~^ eV^. In the case of Solution 



2, we restricted both mass parameters to the same intervals as in Sec. [V, i.e., to 10 '^ eV^ < Am^ < 10 '^ eV^ and 
1 eV^ < AM2 < 10 eV^ 

Figure ^ shows the best point values for the five parameters, which were obtained for the different values of the 
flux uncertainty factor rj in the case of Solution 1. All three mixing angles show the largest deviations from the 
value obtained for 77 = 1 for the smallest values of the flux uncertainty factor, whereas for 77 > 1 they become almost 
stable. Somewhat surprising is the high stability of the small mass squared difference, since it was obtained from 
the minimization procedure with large standard deviations. The large mass squared difference, finally, grows almost 
linearly with the flux uncertainty factor, but remains within a region around AM^ = 10"'^ eV . 

The dependence of the parameters corresponding to Solution 2 on the flux uncertainty factor is shown in Fig. 0. 
The mixing angles are clearly seen to be more stable than in the case of Solution 1, they remain within ranges of 
10° or less around their values at 77 = 1. The small mass squared difference, on the other hand, varies much more 
than the one corresponding to the Solution 1. Finally, the large mass squared difference shows a random distribution, 
as expected, since the minimization function is independent of this parameter in the mass region corresponding to 
Solution 2, as pointed out in Sec. 0. 

IX. SUMMARY AND CONCLUSIONS 

We have fitted the five parameters of a three flavor neutrino oscillation scenario to experimental values obtained 
by the Super-Kamiokande collaboration for atmospheric and solar neutrinos. To obtain numerically stable solutions, 
we divided a large region for the mass squared difference parameters into 45 smaller regions and performed the fit for 
each of these regions. A mass region was considered to contain a solution of the minimization problem, if the best 
point values for mass squared differences were obtained from the fit within the region boundaries. Furthermore, the 
best point value of the minimization function was supposed to go through a minimum in such a region if compared to 
the values obtained in the surrounding regions. In this way, three types of solutions of the minimization problem were 
obtained, out of which the best one corresponds to the most common three fiavor neutrino oscillation scenario, with 
one mass parameter in the "just-so" range and the other one in the interval lO^"' eV^ < AAf^ < lO"'^ eV^. As in this 
common scenario, the first and the third mixing angle of this solution were obtained to be maximal, Oi ~ 45° and 



03 — 45°, respectively, whereas the second one was obtained to be small, 62 — 10°, which is below the CHOOZ upper 
bound. This solution corresponds to the global minimum of the minimization function in the considered parameter 
range, and it was the only one that provides a good fit to all three types of data considered. As a disadvantage the 
value of the large mass squared difference was obtained with a rather large standard deviation. 

The second solution is characterized by a value of the small mass squared difference in the range, which contains 
possible solutions of the atmospheric neutrino problem, i.e., between 10~^ eV^ and 10^^ eV^ but by a large mass 
squared difference in the LSND range, i.e., between 1 eV^ and 10 eV^. Here the mixing angles 6*1 and 82 were 
obtained to be maximal, whereas for the third one there was obtained a smaller value, ^3 ~ 15°. This solution 
provides a comparably worse fit to the atmospheric electron neutrino data, which is the reason for the larger value of 
the minimization function corresponding to that solution. 

The third solution corresponds to a small mass squared difference in the range 10^^ eV^ < Arri^ < 10^^ eV^ 
and a larger one in the LSND range. Here all three mixing angles were obtained quite close to each other, 9i ~ 
38°, 02 — 29°, 03 ~ 25°, the latter one with a rather large standard deviation. This solution provides good fits 
to the atmospheric electron neutrino data and the solar neutrino data. However, the oscillation probabilities for 
the atmospheric muon neutrino data corresponding to this solution show no zenith angle dependence such that this 
solution can be ruled out by the results obtained by Super-Kamiokande. 
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Best point values for AAf^ 
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TABLE I. Best point values of the mass parameters in the various mass regions. All numerical values are given in units of 
eV^. 
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TABLE IL Lines of values for the mass squared differences obtained within the boundaries. A capital "M" denotes a value 
for AM^ obtained within the boundaries, whereas "m" stands for a corresponding value for Am^. The regions labeled by 
"M m" contain solutions of the minimization problem. All numerical values are given in units of eV^. 
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TABLE in. Best point values of the minimization function in the various mass regions. The mass squared differences Am'^ 
and AM^ are given in units of eV^ . 
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Solution 1 
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TABLE IV. Parameter values for the three obtained solutions. 




FIG. 1. The dependence of the ratio Re on cos 6*. 




FIG. 2. The dependence of the ratio R^ on cos 5. 
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FIG. 3. The probability Paun as a function of the solar neutrino energy E. 
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FIG. 4. Stability of Solution 1 with respect to variation of the flux uncertainty factor rj. 
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FIG. 5. Stability of Solution 2 with respect to variation of the fiux uncertainty factor rj. 
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